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Abstract. We consider a consistent linear effective vielbein matter coupling without in¬ 
troducing the Boulware-Deser ghost in ghost-free massive gravity. This is achieved in the 
partially constrained vielbein formulation. We first introduce the formalism and prove the 
absence of ghost at all scales. As next we investigate the cosmological application of this 
coupling in this new formulation. We show that even if the background evolution accords 
with the metric formulation, the perturbations display important different features in the par¬ 
tially constrained vielbein formulation. We study the cosmological perturbations of the two 
branches of solutions separately. The tensor perturbations coincide with those in the metric 
formulation. Concerning the vector and scalar perturbations, the requirement of absence of 
ghost and gradient instabilities yields slightly different allowed parameter space. 
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1 Introduction 

The existence of a graviton mass is an unavoidable fundamental question from a theoretical 
perspective. The pioneering work of Fierz and Pauli marked an important progress [1]. It is 
the unique mass term at the linear level, which does not lead to the presence of ghosts at the 
classical level in the theory. Even though theoretically completely consistent, this theory of 
massive gravity suffers unfortunately from the vDVZ discontinuity [2, 3], i.e.. General Rela¬ 
tivity predictions are not recovered in the limit of massless gravitons. The main reason for 
this is the existence of an additional scalar degree of freedom in massive gravity that couples 
to the trace of the energy momentum tensor. Obviously the vDVZ discontinuity is just an 
artifact of the linear approximation. The effects of massive gravity might be cloaked by the 
Vainshtein mechanism, where the helicity-0 mode interactions become appreciable to freeze 
out the held huctuations on small scales. Nevertheless, such non-linear extensions usually have 
the Boulware-Deser ghost instability when non-trivial backgrounds are considered. Recently, 
a ghost-free non-linear theory of massive gravity was pushed forward [4-7], which has sparked 
a renewed interest in massive gravity, specially for its possible application in cosmology. 

The phenomenological exploration of the theory has triggered a knock-on effect for ex¬ 
tensions of the standard formulation. First of all a theoretically well expected no-go result was 
faced soon [8]. The same condition to remove the Boulware-Deser ghost enforces the absence 
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of flat FLRW solutions in the case of flat fiducial metric. Although a self-accelerating open 
FLRW solution exists [9], it suffers from non-linear ghost instability [10, 11]. Trying to bypass 
this difficulty has initiated considerations of more general hducial metrics, which sadly was 
soon doomed to have Higuchi type instabilities in the case of dS reference metric [12, 13] and 
to protect against acceleration in the case of AdS reference metric [14]. Moreover, the self- 
accelerating branch solutions are not affected by the choice of the reference metric and they all 
behave like the open solutions above [15, 16]. This on the other hand has launched extensions 
of the theory by adding new degrees of freedom [17-19]. The study of quantum properties 
of the potential interactions has motivated a new branch of research for consistent matter 
couplings in the theory. The requirement of maintaining the special structure of the potential 
interactions intact at the quantum level yields severe restrictions on the possible couplings. 
The phenomenological consequences of massive (bi-)gravity depend highly non-trivially on the 
way how matter helds couple to the massive graviton [20-40]. The above mentioned difficulties 
in the cosmological setup arise only in the case of the minimal coupling to one of the metrics, 
which can be avoided by considering a matter coupling through a very specihc admixture of 
the dynamical and hducial metric [25]. Considerations of this direction usually reintroduce 
the Boulware-Deser ghost. The important question there is at what scale this ghostly degree 
of freedom enters. In the context of the special composite effective metric proposed in [25], the 
hrst non-trivial non-vanishing leading order vector-scalar-matter interactions typically rein¬ 
troduce the Boulware-Deser ghost assuming full local Lorentz symmetry [41]. However, this 
does not preclude us from considering these interactions as a consistent effective held theory 
with the cut-oh given by the mass of the ghost. Not only avoids this specihc coupling the 
no-go result for the hat FLRW solutions, it does it so in a very specihc way. The mass of the 
Boulware-Deser ghost on the FLRW background is inhnite. This is because the ghostly vector- 
scalar-matter interactions do not contribute to the background due to the symmetry of FLRW. 

Of course it would be more tempting to construct matter couplings in which the absence 
of the Boulware-Deser ghost is realised fully non-linearly. Not only would this enable the via¬ 
bility of the Vainshtein mechanism, but also enlarge the scale of applicability. A possible way 
out was argued in [42] by switching to the unconstrained vielbein formulation. The argumen¬ 
tation is simple and elegant. In the case of non-derivative coupling the resulting Hamiltonian 
is linear in the effective lapse and shift. In the boosted ADM decomposition of the vielbein 
one can easily express the effective lapse and shift in terms of the lapses and shifts of the two 
vielbeins. After using the equations of motion for the shift one can integrate out the boost 
parameters resulting in effective lapse and shift that depends only linearly on the lapses and 
shifts of the two vielbeins. Hence, the linearity in the lapse and shift enforce the hrst class pri¬ 
mary constraints that remove the Boulware-Deser ghost. So far so good, this argumentation 
has however an unfortunate loophole. The existence of the secondary constraints is taken for 
granted. Nevertheless, this turns out to be not the case. In [43], it was shown explicitly that 
the Hamiltonian becomes highly non-linear in the lapses once the rotations are integrated out 
as well, which rehects the absence of the secondary constraints along the line of the analysis 
performed in [42]. Thus, also in the unconstrained vielbein formulation the Boulware-Deser 
ghost remains persistent in the non-minimal coupling. 
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In this work we will consider yet another formulation of the coupling, which we call 
partially constrained vielbein formulation. For this purpose we shall go along the lines of a 
recent study in [44], where the goal was to remove the unwanted, unstable degrees of freedom 
of dRGT theory, by only keeping its tensor modes. The advantage of this formulation is 
the absence of the Boulware-Deser ghost to all orders beyond the decoupling limit. Any 
general vielbein can be decomposed into a Lorentz boost and rotation of a triangular vielbein. 
In [42], the integration of the {D — 1) boost parameters in D dimensions yields a linear 
Hamiltonian, however the integrations of the additional {D — 1){D — 2)/2 rotation parameters 
yields a Hamiltonian highly non-linear in the lapses [43]. Our partially constrained vielbeins 
are constructed such that the rotation parameters do not reintroduce the non-linearities in 
the lapses. Hence, the argumentation of [42] applies exactly and the loophole found in [43] is 
removed. We shall further study the cosmological consequences of the partially constrained 
formulation of the matter coupling. For this we will assume the same background evolution as 
for the metric perturbations with vanishing boost parameters. However, important differences 
arise at the level of the perturbations. We shall first present the dRGT theory in the partially 
constrained formulation in Section 2 and state our convenient conventions and notations. We 
then prove the absence of the Boulware-Deser ghost related to the non-minimal coupling in 
this new formulation in Section 3. As next we study the background equations of motion on 
FLRW space-time in Section 4 with the presence of two branches of solutions. Similarly as in 
the previous studies of the effective coupling the new coupling enables us to avoid the no-go 
theorem for flat FLRW solutions. Finally, we shall investigate in detail the stability of the 
perturbations in the two branches in Section 5. 

2 Partially constrained vielbein formulation 

The original formulation of de Rham-Gabadadze-Tolley (dRGT) massive gravity in the metric 
language is mathematically cumbersome due to the presence of matrix square roots. Not only 
is the Hamiltonian analysis hindered but also the cosmological perturbations. For instance 
to prove the existence of the primary constraint one has to perform a highly non-trivial held 
redehnition of the shift. Since the vielbein is like the square root of the metric, the study of 
the potential interactions simplify drastically when one uses the vielbein formulation. We can 
express the two metrics by the corresponding vielbeins as 

/.jS u and f^i, tjabE y. (2-1) 

Hereafter, indices A, B, C, V and /r, u, p, a run from 0 to 3, but indices /, J, K, L and i, j, k, I 
run from 1 to 3. In terms of this Einstein-Gartan formulation, the potential interactions 
simply correspond to all possible wedge products of the vielbein of the dynamical metric with 
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the vielbein of the second metric 


r _ ^ tp-A tpB jpC TpV 

'~'0 24*^ ^ABCD-t^ ^h/ i/hj phj fj , 

£2 = \e^''^^ej,scvE-^pE^,e‘^pe^„, 

£3 = ^e^^''P-eABCvE^pe^ue^pe^., 

Ea = pB^cr , ( 2 . 2 ) 

with the Levi-Civita symbol normalized as €0123 = = 1. The Einstein-Hilbert term 

in pnre General Relativity written in the vielbein formnlation is trivial in the sense, that it 
corresponds to introdncing new non-dynamical helds together with gange invariances, exactly 
in the same philosophy as the Stiickelberg trick. It admits fonr diffeomorphisms and six local 
Lorentz transformations. Introdncing the potential interactions breaks these symmetries. In 
bigravity, the two copies of local Lorentz transformations and diffeomorphisms are broken 
down to a single copy. Of conrse one can introdnce Stiickelberg helds in order to restore the 
broken Lorentz and diffeomorphism invariance. The vielbein formnlation of the bi-gravity 
theory is dynamically eqnivalent to the metric formnlation if one imposes the symmetric 
vielbein condition, which is a direct conseqnence of the eqnations of motion for the Lorentz 
Stiickelberg held. 

In this work we will be concentrating on massive gravity, i.e. the case with non-dynamical 
hdncial vielbein together with its dnal basis satisfying 

E^pE^‘' = 5;, E-^pEs^ = 5^. (2.3) 

We also introdnce the dnal basis for the dynamical physical vielbein e"^ as 

= (2.4) 

As already stated, the graviton mass term as well as the Einstein-Hilbert action are invariant 
nnder the overall local Lorentz transformation of the two vielbeins. We can hx the gange free¬ 
dom associated with the boost part of the overall local Lorentz transformation by demanding 
that the hdncial vielbein is of the so called Arnowitt-Deser-Misner (ADM) form: 


^ ^ ^ 

^ \M>^E^^EE)- 


(2.5) 


Here, M, M* and EE are the hdncial lapse, the hdncial shift and the hdncial spatial vielbein 
since the corresponding 4-dimensional hdncial metric becomes 

fp^dx^dx’^ = -M'^dt^ + fij{dE + M^dt){dE + M^dt), fij = 5 ijE\E-^^. (2.6) 
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On the other hand, after setting the hducial vielbein to be of the ADM form (2.5), in general 
the physical vielbein cannot be cast into the ADM form simultaneously. Instead, by using 
the known fact that any vielbein can be written in the so called boosted ADM form, we 
parametrize by the physical lapse N, the physical shift A^*, the physical spatial vielbein 
and the boost parameter bj as 

( JV6' + nV, (s[ + e'i («£ + ^) j ’ 

where b^ = 6^'^bj and 7 = \/l + b^bj. This is not a gauge choice nor a physical condition 
but simply a reprametrization of a general vielbein. In terms of these variables the physical 
metric is written as 

= —N'^dt^ + gij{dx'' + N'^dt){dx^ + N^dt), gij = 5ije\e'^y (2.8) 

After hxing the boost part of the overall local Lorentz transformation as (2.5), we further 
impose the following symmetric condition on the 3-dimensional spatial vielbein: 

Yu = Yjj, (2.9) 

where Yjj = E/djxe^i- The condition (2.9) is not a gauge condition but is a physical condition 
since it is invariant under the spatial rotational part of the overall local Lorentz transformation. 
Imposing (2.9) thus modihes the original formulation and dehnes a different formulation. Since 
the physical condition (2.9) treats the spatial components and the temporal component of the 
vielbein in a different way, this formulation potentially violates local Lorentz invariance in 
the gravity sector. We call massive gravity with the symmetric condition (2.9) a partially 
constrained vielbein formulation of massive gravity. 

We thus have (at least) four formulations of the dRGT massive gravity: the metric 
formulation, the constrained vielbein formulation, the unconstrained vielbein formulation and 
the partially constrained vielbein formulation. All of the four formulations of the dRGT 
theory are equivalent on shell, i.e. after imposing equations of motion. This is because in the 
unconstrained or partially constrained vielbein formulation, the equations of motion set the 
4-dimensional vielbein to be symmetric on shell. As we shall see in the following sections, this 
is no longer the case if we introduce composite coupling to matter helds. The advantage of the 
partially constrained vielbein formulation is that, unlike other formulations, matter helds can 
consistently couple to a composite vielbein without turning on a (would-be) Boulware-Deser 
ghost at all scales. 

We have already hxed the boost part of the overall Lorentz transformation by demanding 
the ADM form as in (2.5). On the other hand, we have not hxed the spatial rotational part 
of the overall Lorentz transformation. While it is not necessary, for practical purposes it is 
sometimes convenient to hx the spatial rotational part as well, e.g. by introducing a hxed 
non-singular 3x3 matrix and imposing the condition 

= E\£^^. ( 2 . 10 ) 


5 



Because of the standard polar decomposition (applied to the 3x3 matrix different 

choices of 8'^^ simply correspond to different choices of the spatial rotational part of the overall 
local Lorentz transformation, on which the graviton mass term and the Einstein-Hilbert action 
(as well as the matter action, even with the composite coupling discussed below) do not 
depend. For example, for the analysis of perturbations around a particular background, it is 
probably most convenient to choose 8'’^ to be the background value of 5’^^Ej^. 

On introducing the variables and as 

Ye^ = E/e^^, E-*^, ( 2 . 11 ) 

we can rewrite the potential interactions (2.2) simply as 

Co = -^, 

C-i = 

A = - Y/Y^^) , 

C3 = -V^X^-^, 

A = (2.12) 

In the absence of the composite matter coupling introduced in (2.13) below, i.e. if matter helds 
minimally couple to the physical metric then the symmetric vielbein (Y-^® = Y®-^) is a 
solution to the equations of motion (as in the unconstrained vielbein formulation) and thus 
the partially constrained vielbein formulation agrees with the symmetric vielbein formulation 
and the metric formulation. 

Concerning the coupling to the matter helds, we will only allow non-derivative couplings 
through a specihc composite effective metric proposed in [25, 28]. In terms of the composite 
effective vielbein 

= ae-^^ + (2.13) 

it is constructed as 

— VAB&tstJ.^eSv ) (2-14) 

with constants a and 13. Hence, the action that we will be considering in this work is 

dA [^^ det e R[e] + Mp^rri^ ^ A A + det Ces A(eeff, 0)] , (2.15) 

n 

where 0 represents the matter held. With the composite matter coupling, the partially con¬ 
strained vielbein formulation results in physical consequences diherent from those in other 
formulations, as we shall see in the next section. For the specihc homogeneous and isotropic 
background of interest we will be concentrating on solutions in which Y-^® is symmetric. 
Nonetheless, perturbations around it introduce the antisymmetric part of Y-^^ in diherent 
ways in diherent formulations. 
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3 Absence of the Boulware-Deser ghost 


The Einstein-Hilbert kinetic term for the physical metric is invariant under two separate local 
Lorentz transformations, 

^ ^ (3.1) 

However, the inclusion of the graviton mass terms breaks the two independent copies of local 
Lorentz transformations to a diagonal combination. Hence the action depends on A“^A, while 
it is independent of the overall local Lorentz transformation. Without loss of generality it is 
thus possible to make use of the overall Lorentz transformation to set the fiducial vielbein 
to the ADM form (2.5). Moreover, any general vielbein can be written in the boosted ADM 
form. We thus reparametrize the physical vielbein as (2.7). Similarly, we reparametrize 
the effective vielbein by the lapse A^^eg, the shift A(fg, the spatial vielbein and the 
boost parameter as 



NeSleS + eefffc 




k 




5i + 




7eft+l y effj 


5i+ ’ 


7eft + l ) 


(3.2) 


where = and 7 eff = + b^f^bf^. In the following we will closely follow the argument 

presented in [42] to show that AA.fr and A"®g are linear in N, N\ M and M\ Substitution of 
the above decompositions of the two vielbeins e^, and the effective vielbein into the 
relation + f3E-^^ results in 

abie^j, 

T (3-3) 

and 


UeS^L 




■'eSj 


Sl + 


7efr + 1 


A^effbefif + N^sb'h^eak 

AT hi I Ajk L ( 

JVeffOefr + JVggegg^ I 0^ H-— 

V 7efr + 1 


a[E-i + E^bLe\)+l3M, 

= a Nb^ + N^e\{5i + ^^ 
V 7 + 1 


+ (5M'^E\. (3.4) 


Since the right hand sides of (3.3) do not depend on N, N^, M and M*, by solving these 
equations with respect to and b'f, gives rise to expressions of e^g^ and bf^ in terms of 
e^p, E^p and 6^. By substituting these solutions to (3.4), one obtains a set of linear equations 
for AAff and A^^^g. The right hand sides of these linear equations are still linear in N, N^, M 
and ML Therefore, by solving them with respect to A"efr and A7g, one sees that A"efr and A^gg 
are linear in N, N^, M and ML Furthermore, the Hamiltonian density of the matter fields 
propagating on is linear in A^.^ and A7g. Thus, it is linear in N, N\ M and ML Needless 
to say, the Hamiltonian density of the matter helds propagating on is linear in N and A^* 
and independent of M and ML Also, each term in the graviton mass terms (2.2) includes 
only one factor of or E^q. Since e\ and E^q are linear in N, N\ M and M* and other 
components of vielbeins are independent of N, A^*, M and M®, the graviton mass terms are 
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linear in N, TV®, M and M®. Therefore, the Hamiltonian of the system is linear in TV, iV®, M 
and M®. Following [45], one can thus show that there is a primary constraint that removes the 
Boulware-Deser ghost as follows: The equations of motion for {i = 1,2,3) do not depend 
on TV, M and but depend on the boost parameters hi {L = 1,2,3). Thus one can solve 
these 3 equations with respect to the 3 parameters and the solutions are independent of 
TV, TV®, M and TIV®. Therefore, after substituting them, the Hamiltonian is independent of TV* 
and linear in TV, M and TIT®. The equation of motion for TV is thus a constraint. This is the 
primary constraint that removes the Boulware-Deser ghost, with the help of the associated 
secondary constraint. 

In the case studied in [42], this argumentation unfortunately had a loophole. This loop¬ 
hole was due to the fact that the structure of the Hamiltonian drastically changes after solving 
the equations of motion for the spatial rotational part of the local Lorentz transformation. 
This fact was not appreciated in [42] but was later noticed in [43]. In the metric formulation (or 
equivalently the constrained vielbein formulation) of the dRGT theory without the composite 
matter coupling, the (would-be) Boulware-Deser ghost is removed by a primary constraint 
and the associated secondary constraint in the phase space spanned by gij and their canonical 
momenta. The presence of such a primary constraint is equivalent to the fact that the Hamil¬ 
tonian density is linear in the lapse after integrating out the shift vector. The unconstrained 
vielbein formulation (after hxing the overall local Lorentz transformation), on the other hand, 
includes not only the lapse and the shift but also the boost parameters and the spatial 
rotation parameters as basic variables. One thus needs to show that the Hamiltonian density 
is linear in the lapse after integrating out not only the shift but also the boost and spatial 
rotation parameters. This is indeed the case in dRGT theory without the composite matter 
coupling. In the presence of the composite matter coupling, the argument [42] proves that the 
Hamiltonian density is linear in the lapse after integrating out the shift and the boost param¬ 
eters. However, there still remain the spatial rotation parameters. However, as pointed out 
in [43], the Hamiltonian density becomes highly non-linear in the lapse after integrating out 
the spatial rotation parameters as well. Thus the unconstrained vielbein formulation of the 
effective coupling reintroduces the Boulware-Deser ghost in analogy to the metric formulation. 
Gonversely, in the case of partially constrained vielbein formulation that we study here, the 
argumentation of [42] is solid. The spatial rotational part of the local Lorentz transformation 
is already hxed by the symmetric condition (2.9) at the level of the theory and thus there 
are no spatial rotation parameters in the set of basic variables. After using the equations of 
motion for the shift, we can integrate out the shift and the boost parameters all together, and 
the Hamiltonian density remains linear in the lapse as explained above. Hence the loophole of 
the argumentation in [42] is removed by the symmetric condition (2.9), at the price of losing 
local Lorentz symmetry. Our setup illustrates a correct realization of the argument given in 
[42]. 

4 Cosmological background evolution 

The no-go result for flat FLRW solutions in the standard formulation of massive gravity [8] has 
motivated to consider non-minimal matter couplings. In the standard coupling the Stuckelberg 


field equation of motion imposes a non-dynamical scale factor. Non-minimal matter coupling 
through a consistent composite effective metric softens the restrictions on the scale factor 
and one can easily construct exact FLRW solutions with flat reference metric [25]. In the 
vielbein formulation this effective composite metric corresponds to a linear effective vielbein 
built additively out of the two vielbeins. In the metric formulation this non-minimal coupling 
reintroduces the Boulware-Deser ghost, however, this does not prevent one to view this cou¬ 
pling as effective field theories below some cutoff scale. It was argued that the Boulware-Deser 
ghost might remain absent in the unconstrained vielbein formulation, but unfortunately the 
absence of the secondary constraints was proven very soon. The setup that we are using here 
corresponds to the correct realisation of the argument given in [42]. In the previous section 
we have explicitly shown the ghost absence in the partially constrained vielbein formulation. 
This motivates us to investigate the cosmological implications of this model further. For this 
purpose, we shall assume a background with = 0. Note that from the graviton mass 
terms one only obtains non-linear contributions in Thus, in the absence of the matter 

coupling through the linear effective vielbein one consistent solution would be F[^b] = 0. Even 
in the presence of the matter coupling we shall concentrate on the background with symmet¬ 
ric Y_aj 3 for the FLRW background. Let us consider the following homogeneous and isotropic 
co-diagonal backgrounds for the dynamical and hducial metric. In other words both metrics 
have a FLRW form in the same coordinate system ^ 


ds^g = —N'^dt^ -|- a^dijdx^dx^ , 

ds^f = —N'jdt^ + ajdijdx^dx^ . (4.1) 


The corresponding tetrads read 


e 


0 




K = > 



and similarly the diagonal symmetric matrix Y] 4 b takes the form 


Yab = diag 


f N a a a \ 
y N f Ojj Cij- 


(4.2) 


(4.3) 


Accordingly, the composite effective metric on this background corresponds to the line element 

dSgg = -Nl^dt^ + algdijdx'dx^ , (4.4) 

where Nes is the effective lapse function and aes stands for the effective scale factor 


Nes = (y N + (3 Nf , ttes = C(a + (3 ttf . (4.5) 

^ Although we use a generic FLRW form for the fiducial metric, the geometry of the fiducial space in 
massive gravity should be specified at a fundamental level. This connection with a fundamental geometry 
might then impose further restrictions on Nf and a/. For a list of FLRW charts of maximally symmetric 
space-times, see e.g. Ref. [16]. 
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As for the matter sector, we will only consider a generic scalar field 0 minimally coupled to 
the effective metric 


S<j) 


j d'^xdet Ceff £,/,(eeff,0) 


d'^x det CefT P (x, 0) , 


(4.6) 


with X = — Since the scalar held has to be compatible with our above homo¬ 

geneous and isotropic Ansaetze, we will impose that the scalar held depends only on time 
0 = 0(t) at the background level. For later convenience let us also introduce the following 
quantities 


P= “ ,Hf- , 

a AP ^ a/ iV/ 

(4.7) 

_ Nf/af 

N/a ' 

(4.8) 

U (A) = 0i4A^ -|- 602A^ + 403A, 

(4.9) 


with H denoting the expansion rate, A = oo/a the ratio of the scale factors and r the speed of 
light propagating in the hducial metric with respect to the one propagating in the dynamical 
g metric. The action on the homogeneous and isotropic background simplihes to 


^ y dt I - A - - 771^04 - [pm + rAQ] | 

J dt(3oajNf + J dt al^NesP {x, (p) , (4-10) 

where we introduced the shortcut notations for manageability 

P^(A)^U{A)-^Ua, Q{A)^^Ua, (4.11) 

with U^A = SaU. After varying the mini-superspace action (4.10) with respect to the lapse N 
we obtain the Friedmann equation 

= A + + rn^pm + jt^P , (4-12) 

iVj pJ CL 

where it becomes clear that pm denotes a dimensionless ehective energy density from the mass 
term and p = 2 0^P(x, 0)x “ P{Xy<P) is the associated energy density of the matter held. 
The acceleration equation results from combining the equation of motion for the scale factor 
a with the Friedmann equation 


2H 


J A{r — 1) 





clNqs p 


(4.13) 
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with J{A) = ^dAPm{A) and the pressure of the matter held P = P(x,0)- Similarly, the 
variation of the mini-superspace action with respect to the scalar held yields 


1 


p + 3 


^efF 

Oefr-^efT 


(p + ^) = o, 


(4.14) 


which just correspond to the standard conservation equation for a held minimally coupled to 
the Peff metric. Finally, the contracted Bianchi identity gives the constraint equation 


{ H - H,A) 


MpiJ 


a 13 als 


0 . 


(4.15) 


Generically, this equation dehnes two branches of solutions: Branch-I H = Hf A and Branch- 
II J oc P. When we study the perturbations, we will consider both branches of solutions 
separately. At this stage, it is worth to mention that for Minkowski hducial metric a/ = oq 
and Hf = 0, only the second branch gives an expanding solution. In fact integrating the above 
equation for the case of Minkowski hducial metric yields [46] 

Q + (3algp = alm^Mp^K , (4.16) 


where k is a dimensionless integration constant, independent of the normalization of the scale 
factor. The contracted Bianchi identity together with its integrated version in the case of 
Minkowski hducial metric permit to relate the kinetic and potential energy density of the 
matter held to the graviton mass terms. The Friedmann equation can thus be brought into 
the form 


3P2 


A -I- + m? 





(4.17) 


Similarly, one can also remove the explicit y dependence in the acceleration equation (4.13), 
yielding 


2H 


= —m 


JA-^{Q-kA^-J) 


(4.18) 


which is nothing else but the time derivative of the previous equation (4.17). In the present 
study, we will not make use of this fact, and we will keep the hducial metric arbitrary. 

As expected, the background evolution is exactly the same as in the metric formulation. 
Note that this is not surprising since we considered here background solutions with symmetric 
Yab- There could be other interesting background evolution with a non-zero YyAB] background 
values. Even in our case with the vanishing Y[ab] we expect non-trivial diherences to the metric 
formulation at the level of the perturbations. 

For later convenience in the study of perturbations we introduce the quantity 


T = AJ + 


(r — 1)A^ 
2 


OaJ, 


(4.19) 
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and the following relations that hold 


—pm = 3JA{HfrA-H), 


1 • 2{HfrA 
N'^ ~ A{r-1) 

L = ^HfrA-H) 


H) 


{T-JA) , 


(r - l)pm + ^ ^ r 


3 — 2 r + 1 


JA 




5 Stability of the perturbations 


As we have seen in the previous section, the coupling of the matter helds through a very 
specihc admixture of the dynamical and hducial vielbeins avoids the no-go result for flat 
FLRW solutions in massive gravity. In the metric formulation it was explicitly shown, that 
the Boulware-Deser ghost is absent around FLRW solutions [25]. In our partially constrained 
vielbein formulation this is true in general for any background. The ghost is maintained absent 
non-linearly by breaking explicitly the local Lorentz invariance. Even if the Boulware-Deser 
ghost is removed, it is important to investigate how the remaining physical degrees of freedom 
behave on top of the background that we are interested in. For this purpose we shall adopt 
the boosted ADM formulation where the physical tetrad can be written as 


where 


ek = (e-“)k 

are the Lorentz transformations and 


B 

At ; 

corresponds to the ADM tetrad 


(5.1) 


(5.2) 

with the lapse perturbation <F. For the partially constrained tetrad formalism, only the boost 
part of the Lorentz transformations are kept, i.e. uu = 0. The perturbations of the boost 
parameters are decomposed into 

cjQi = div -f Vi , (5.3) 

where 6^'^djVi = 0. The induced 3-metric can be built out of the triads eh as hij = ^dij. 
The shift vector is decomposed as 


e-^„ = 


iV(l + <F) 0 




= = aN{Bk + dkB), 


(5.4) 


while the triads are decomposed through 


= a (1 + 'ip)5l + 




lij + d{iEj) + 


5^6^- 


-5. 
3 


5 "' 1 d,diE 


V'- 


(5.5) 


where 5^WiBj = S^^diEj = = 0. Notice that the decomposition of the metric 

perturbations coincide with the ones in [46] at linear order. Finally, the perturbation to the 
scalar held is introduced by 


(j) = (j)o + S(j). 


( 5 . 6 ) 











In the following we stndy in detail the stability of tensor, vector and scalar perturbations in the 
two branches of solutions and exploit the conditions to avoid ghost and Laplace instabilities. 
For this purpose we will further decompose the perturbations in Fourier modes with respect 
to the spatial coordinates 

/ d^k 

exp{ik ■ x) + c.c., (5.7) 

where T represents the perturbations ($, E, B, 6(j), v) and (i?/, Ej, vj, hjj) respectively. 
In the remainder of the paper, we will omit the subscript k associated with the mode functions 
for the sake of clarity of notation. We shall perform the stability analysis of the perturbations 
for each sector separately in what follows. 


5.1 Branch-I 

As we mentioned above from the constraint equation (4.15) one obtains two branches of 
solutions. Let us hrst investigate the perturbations on top of the Branch-I solutions with 
H = AHf and analyse the parameter space that is free from the ghost and Laplace instabilities. 

5.1.1 Tensor perturbations 

We start the stability analysis with the transverse-traceless fluctuations. Our action (2.15) 
perturbed to second order in tensor perturbations has the form 


(2) _ M|i 

- 


■^tensor 


d'^k dt N 


— - 
N2 IJ 






with the mass term dehned as 


2 _ 2-n 

rrirp = m L 


a P ttesNesP A 
Mp^a N 


(5.8) 


(5.9) 


First of all the tensor perturbations have the right sign for the kinetic term and also the right 
condition for absence of gradient instabilities. The only condition comes from the requirement 
to avoid tachyonic instability rriT > 0. 


5.1.2 Vector perturbations 

As next we can compute our action perturbed to second order in vector perturbations: 


5, 


( 2 ) 


vector 


= —[ d^k N dt k^a^ 
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aN 


E*jB^ + BIE^ 


mlE*jE^ 




ml + — mis {v]B^ + B*v^) 


' k'^ 

where we have introduced the following dehnitions for convenience 


(5.10) 


ml =m‘^JA{r + l) + 


a jd A a \ 


2 

eff 


M2ja2 


(r+l)p-L^(p + F) 

aeffA 


m. 


2 _ 
vB 


= m J A -|- 


a(d A OgffP 


(5.11) 
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Note that the vector helds vi and Bj are non-dynamical degrees of freedom. We can therefore 
compute their equations of motion and integrate them out. Solving the equations of motion 
for Vi and Bj, we obtain 


m. 


Vi = 


vB 


mt 


B, 


Bj = 


2N 


1 + 


2 a? mt 


vB 


A;2 


mt 


-1 


Ej 


We can plug back these solutions in the quadratic action (5.10) and obtain 


5, 


( 2 ) 


vector 


Mi 

16 


d k N dt k a 


2„3 




EJE' - mlEJE’ 


(5.12) 


(5.13) 


1 + I ■ have to impose that the vector perturbations have the right 

sign for the kinetic term /Cy > 0 in order to avoid ghost instabilities. We can rewrite the 
kinetic term explicitly as 



1 + 


2Aa^ 


r + 1 


^2 J I 

m J + ^2^ 


® eflf -^eff 

oivM 


a P {p + P) 


m'^J + 



(5.14) 


At high momenta, the condition is rather cumbersome. At low momenta in the regime where 
MppnPj 3> p, P, the kinetic term is positive if m? J > 0. In the other regime, where the 
matter is parametrically dominant over the mass terms, i.e. MpprAj p, we get the 
condition ^ 

P > ..rffr+il ■ (5-15) 

a VefI 


5.1.3 Scalar perturbations 


As next we will study the stability of the scalar perturbations in our partially constrained 
vielbein formulation of the dRGT theory with the effective coupling. As we mentioned above, 
hve dof appear in form of a scalar held <I>, B, v, 'ipi E, 6(j), from which three (<I>, B and v) 
are non-dynamical and will be integrated out. The equations of motion for $, B and v are, 
respectively 


a N 

+-Ai’ + -E 


1 

TV 


+ 


al^N (p + P ) ^ , 3A 

2 M2ja3iVefrc2 
aaeff(P + ^) 


<h + M + 





2 M|[ PclAo 

(^alaiP + P) 


5(p + 


a a 


eff 


2 M2 


P 


A 


Pit 


a A 


2 M2j a2 


m^J + 


a A QeffP 

a2M2j 


B - 


(r -I- 1) m J + 


(j>o/NeS 
a fd alup 


— \mJ + 



a?Ml^ 


(50 = 0, 


n = 0, 


a2M2j 


Q(/?aeffAeff(p + P) 


M2jaiV 


ajdaenir - l){p + P) 
M|ia2(^o/^eff 


(50 = 0, (5.16) 
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where we made use of the sound speed of the matter held 


2 _ P,x _ 4>) 

~ p.x '28lP{x,4‘) X + dxP(x,4>)' 


(5.17) 


Using the solutions of these equations reduces the action to contain three degrees of freedom 
only. We hnd that a convenient basis to remove the would-be Boulware-Deser degree is 


Id 


(50- 


a (do 
hfiVeff 




(5.18) 


This choice of basis is just convenience since the subhorizon limit of the kinetic term turns 
out to be of order k^. The basis used for Branch-II will be a different one. In this basis the 
mode 0 becomes non-dynamical. Integrating it out, we end up with two dynamical modes, 
with action 


( 2 ) ^ ^P 1 

^scalar o 


(fkNdt a? 


dv 


y yt y 

■K- — + — -M-Y-Y^-M - Y^ -Xf -Y 

N N N 


(5.19) 


The exact form of the 2x2 matrices = K, = 12^ and = —M are very 

cumbersome, although the subhorizon expansion yields relatively manageable expressions. 
The components of the kinetic matrix at leading order in this expansion {0{k^)) gives 


Kn 

K 22 


Ku 


alsNN,s{p + P) 
00 


+ 0{k 


-ii 






2 iJ 2 


r (J -|- 


a 2 M^j 


a^Mpj J 


a^al^NAip + P) 




+ 0{k 


- 1 ) 


g algN A{p + P) 
2Ml^aclH4)Q 


m^J + 




0{k 


-ii 


(5.20) 


One can rotate the basis such that the eigenvalues become 


Kl — it'll 


alsNN,s{p+P) 

(dg 


+ o{k-^), 


^2 


det K Sa'^A f ^2 j 1 ^ '^esP\‘^ 
Ku ^ 2r(3 ^ a?Ml^ ) 


Q-eff 








“/5“eft(P+'P) j 


— a 




.(5.21) 


In the regime where the matter sector is parametrically subdominant under the mass term, the 
second eigenvalue gives a condition reminiscent of Higuchi bound, i.e. rn^J{m?J A — 2H^)/r > 
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0. In massive gravity, Hf itself is determined by the geometry of the hducial metric, which 
is a fnndamental quantity. Therefore, keeping Hf generic implicitly neglects some necessary 
information. For instance, if the hducial metric is maximally symmetric, this gives a new 
relation between J and the matter components. For this example, the second kinetic eigenvalue 
is simply 


3 a + F) (2 - 1 ) + ' 

2 M^i(r - 1)2 (2 - 1) + 


(5.22) 


which is positive if the sufficient conditions a > 0 and r > 1 are satished. This example can 
also be applied to the vector mode and the vector kinetic term can be found to be positive if 
these conditions, along with f3 > 0 are satished. The anti-symmetric mixing-matrix is of order 
Mi 2 = 0{k^), so it does not contribute to the speed of propagation. Finally, the components 
of the mass matrix are 


(f^^)22 


aeS N^sip + P) 

a3M2jiV02 


+ 0{k ), 


—2 -|- 


A + 


J 


2 



aalgNes{p + P) 



2 I 2 



o 2 \ 2 

a0a^\ 

a^Mpj J 


4 


/ a^aesNesip + P) 

V M2jaih2iv 



-|- 0{k ), 


_ aa.,N^A{p + P) ^ 


(5.23) 


For a mode with subhorizon frequency ui = Cs - + 0{k^), and taking into account the leading 
order expressions of the above matrices, we can obtain the propagation speeds of the modes 
Cg by solving 

det [— K + = 0 , (5.24) 

or 4 ^ 

^ det[JF] (Tr[f22 . R] - Trifl^lTriJF]) Cl + detifl^l = 0 . (5.25) 

o'* H 

5.2 Branch-II 

We would like now study also the behaviour of perturbations on top of the Branch-II solutions 
with the constraint equation J = Similarly as in the previous section we will 

TTX IVjpjCt 

perform the stability analysis of the perturbations for each sector separately. 
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5.2.1 Tensor perturbations 

We shall again start our stability analysis with the transverse-traceless part of fluctuations. 
The action (2.15) perturbed to second order in tensor perturbations is this time 


( 2 ) ^ 
^tensor o 


d^k dt N 


N 2 IJ 





(5.26) 


with the different mass term dehned now as 


O jd Ueffeff P 

) ■ 


(5.27) 


Again the tensor perturbations on top of the Branch-ll solutions have already the right sign for 
the kinetic term and also for the absence of gradient instabilities. The condition for absence of 
tachyonic instability this time becomes fhrp > 0. Actually, the tensor perturbations coincide 
exactly with the tensor perturbations obtained in [46] for the metric formulation. This is no 
surprise, since the boost parameters cuyig contribute only to the vector and scalar perturbations 
(5.3) in our partially constrained formulation of the vielbein. Even if the mass term has 
a very different form at hrst sight, it is actually the same as in Ref. [46]. This becomes clear 
after using the following identihcation of the parameters between the two formulations: 

3 3 

fdi —y —2 ^ ^ -l- 3 0:3 -|- 6 0:4 , —y ——(2 02 + 3 0:3 -I- 4 0:4) . (5.28) 


With this identihcation, we observe that J does not get modihed from one construction to 
the other, while 


(4-1)3 p-l)"[(2r-l)4 + r-2|J-(r-l)(Q + 4y„)] . (5.29) 


Using these relations, equation (5.27) reduces in fact to the one in Ref. [46]. 


5.2.2 Vector perturbations 

As in the previous section our vector modes consist of the non-dynamical degrees Bj, vj and 
the dynamical degrees Ej. The action quadratic in vector modes reads: 


‘^v^tor = I d^k N dtk'^a^ 


16 


- — ( E*B^ + B*E^ ] - fnlE*E 


N^- 


aN 


* TdI I D* ipl \ ^2 7TI* Tpl 


LT-L^I-i 


+ - --ml vW + - mlj, {v]B^ + B^) 


k^ 


(5.30) 


where this time 


_2 a/3 Aaes{ar +/3 A) {p + P) 

'^1, = -- 

a Mpj 


_ 2 _ a^Aalff{p + P) 

^vB — ~ 




(5.31) 
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We remark that the metric formulation corresponds to u/ = + r), reducing to the action 

in [46] at this stage. This can be verihed by looking at the mass terms £„ quadratic in 
perturbations. In the boosted ADM formulation, one has 




Svi 


= 0 . 


vi=Ni/{l+r) 


(5.32) 


In other words, the fully constrained tetrad formalism, which is equivalent to the metric 
formulation, requires ujij = 0 and ujqi = Ni/{1 + r) at linear order around FLRW. This is 
nothing but the symmetric vielbein condition, i.e. Yiab] = 0- We again depart from the metric 
formulation by solving the equations of motion for u/, together with Bj, obtaining 


m. 


Vi = 


vB 


Bi, 


Bj = 


2 \ -1 


1 + 


2 a‘^ rri. 


vB 




Er. 


Using these solutions back in the action, we obtain: 


42 tor = ^ / d^kNdtk^a^ 
16 ./ 




The kinetic term for the vector mode can be written explicitly as 


[I + 


k"^ ml ^ ^ 


2a^mtB 


= 1 + 


e Ml^{ar ^ 13 A) 


-1 


(5.33) 


(5.34) 


(5.35) 


2 a j3a? A{a3r (3 A)^ (P + p) 

At high momenta, assuming p + P > 0, the no-ghost condition corresponds simply to a > 
0,/3 > 0. 


5.2.3 Scalar Perturbations 


As we introduced in the previous section the scalar sector contains six degrees of freedom in 
total: <h, B, ijj, E, v and 6(j). Out of these <|), B, v and a linear combination of the remaining 
ones are non-dynamical. Again, the metric formulation corresponds to hxing the boosts with 
V = B/{1 + r). The action, as usual, is not suitable for presentation. However, we quote here 
some of the intermediate steps. The equations of motion for $, B and v are, respectively 


H 


1 

N 


a N 


+ 


aalfi{p + P) faN^ 3f3aA'ip 


2 


Co A^eff ^eS 


a A'eff 
2 M2 a3 00 




0 + —E - 


2a{H- HfV A) 
cl alsir -1)A 


aals{p + P) 

2 Mpj a3c20o 


Oi [3 OefrA^efF 


(50 = 0, 




+ 


aal^{p + P) { (50 

2Mpja2 V0o/Aefr 


OefT B — a{ar -\- (3 A) V — = 0 . 

0o/Aefr 


+ 3 Hes(p + P) — 


p. 




dcf) 


(3 aAv] = 0 , 


(5.36) 
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Using the solutions of these equations reduces the action to contain three degrees of freedom 
only. The convenient basis to remove the would-be Boulware-Deser ghost on top of Branch-II 
solutions is this time 


Ui = <50 - 


a (j)Q 
HK 


eff 


0 + -U 

D 


Y, = -rE. 


(5.37) 


Once 50 and E are expressed in terms of Yi and 50, the mode 0 becomes non-dynamicah 
Integrating it out, we end up with two dynamical modes, with action 


^ / df’kNdta^ 


yt _ y yt _ _ Y 

—-K- —+ —-M-Y-YYm - yt-o^.y 

N N N N 


(5.38) 


with 2x2 matrices = K, and = —M. 

k‘^{ar + (3 A) — 3 (3 A ^ 


Ki2 = 


Ku = 


2 k NeS H 

a <^o 

KuH 


-1 


k‘^{a r + 13 A) —313 A^ 


Koo = 


kal3a?A^^ 
Ku 


K 


11 


Ki2 + 


Veff N 

2kH 


a 


<t>Q 


k(3a^A%(2kH + K^2a^) 


H 


k‘^{ar + (3 A) — 3l3a?A^ 


(5.39) 


where we have dehned 


H 


3 a 13 cl{H - HfA){r - 1) A 
Ogg {(3 H + a Hf) rfi^ 

One can rotate the basis such that the eigenvalues are: 


-1 


(5.40) 


Ki — Kii — 

det K 

K,2 = 


-{ar + (3 A) — 3 (3 A 


a^/3A 


1 - 


k‘^{ar + (3 A) — 3 (3 a? A 


-1 


K 


= 2E 


11 


H \ ( <j>o\ 

NH^J \N,b ) 

— ^- {ar + P A) 

a3HI3A \ 


a {3 a^A 




^<3 




(5.41) 


The full no-ghost conditions are quite opaque, so we consider the subhorizon limit of the 
action. In this limit the kinetic matrix is diagonal, with K 12 = 0{k~^) and 


i^ii = 


K 22 = 


al^{p + P) 

IVlpU Cg 

a (3 al^A {p -h P) 
Mpja(a r + (3 A) 


1 + 


3 (3^CgiX — r) A H — HjA p-|-P 


2 a 


l3H + aHf 


+ O(fc-0, 


+ o{k-^). 


(5.42) 
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The mixing matrix already has one independent component, which is, at leading order: 


Mi2 — 


aDalg{Hr-H,A)(p + P) 

H + a Hj)(ar + 13 A)-^ 


k + Oik’- 


(5.43) 


Finally, the mass matrix is also diagonal at leading order, with 

,25 (P + P) ,2 




a^M‘^{ar + (3 A)^ 


-e + 0{k), 


(fi')22 = 


2m^{H - Hf A) Nes 

3N A{r-1){^H + aHf) 


e + 0{k). 


(5.44) 


Then, for a mode with subhorizon frequency uj'^ = + 0{k), we can obtain the 

propagation speeds of the modes C| by solving 


det 


uj'^K+ { —K + 3HK + 2M] {-iu) + [SHM + —M + n 


N 


= 0 . 


(5.45) 


Considering the leading order terms of the matrices, the equation can be simplihed as 

det [—K + 2 M{—i ui) + = 0 , (5.46) 


— KuK 22 Cl- — (4 + iF22fin) Cl + = 0 (5.47) 

Qj Qj 

This is a quadratic equation for C| and the solutions can be easily found, although they are not 
very instructive. An interesting limit is /3 —?■ 0, which corresponds to the minimally-coupled 
matter limit; 


CL = cl + 0{f3) , 


^2 

'^S,2 


2 M2j r {H 


HfA) 


m% 


3a'^Hf{r-l)A^{p + P)I3 


+ 0(/3° 


(5.48) 


In other words, for vanishing f3 one of the modes becomes instantaneous, while the other one 
is the matter degree as it propagates with the sound speed of the fluid. This is the standard 
“vanishing kinetic terms” issue of dRGT in self-accelerating branch with minimally coupled 
matter [15]. 


6 Conclusions 

The question of consistent matter couplings within the framework of massive gravity is crucial 
for maintaining the ghost freedom. The specihc quantum properties of the theory have moti¬ 
vated the consideration of an effective matter coupling [25, 47]. This matter coupling through 
a very specihc admixture of the dynamical and hducial metric avoids further the No-go result 
for hat FLRW solutions. However, this coupling reintroduces the Boulware-Deser ghost. It is 
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an indispensable question to address at what scale this ghostly degree of freedom enters. Of 
course these interactions can be still considered as a consistent effective held theory with the 
cut-off scale dictated by the mass of the ghost. On top of a FLRW background the correspond¬ 
ing Boulware-Deser ghost remains absent at second order of perturbations at the level of the 
action and its mass on this background is inhnite. A crucial step would be to study the per¬ 
turbations on top of an anisotropic background and investigate in detail the mass of the ghost. 

Even if one can use the matter coupling through the composite effective metric in the 
sense of effective held theory, it would be more appealing to hnd consistent matter coupling 
in which the absence of the Boulware-Deser ghost is realised fully non-linearly. It was argued 
that this might indeed be possible in the unconstrained vielbein formulation of the theory 
[42]. Using the boosted ADM decomposition of the vielbein one can easily show that the 
ehective lapse and shift functions remain linear in the lapses and shifts of the two vielbeins 
after integrating out the boost parameters resulting in hrst class primary constraints. With 
the help of the associated secondary constraints the Boulware-Deser ghost can be eliminated. 
Unfortunately, the absence of the secondary constraints was shown very soon [43]. In other 
words the integration of the rotations gives rise to a highly non-linear dependence on the 
lapses. This means that also in the unconstrained vielbein formulation it is not possible to 
avoid consistently the Boulware-Deser ghost. 

In this present paper we exhausted yet another formulation of the theory, the so called 
partially constrained vielbein formulation (along the lines of [44]), with the attempt to remove 
the Boulware-Deser ghost fully non-linearly. This formulation guarantees the absence of the 
ghost to all orders beyond the decoupling limit. For this purpose we also adapted to the 
boosted ADM decomposition. After proving the ghost absence, we studied the cosmological 
application of this formulation within the context of the non-minimal matter coupling in 
form of a k-essence model. Even if the background evolution is the same as in the metric 
formulation for vanishing background boost parameters, the dynamics of the perturbations is 
crucially different. This allowed us to put different constraints on the parameters of the theory 
coming from the requirement of ghost and Laplacian instabilities absence. An important 
question to further investigate would be whether or not the ghost absence can be maintained 
in a different decomposition of the vielbein rather than the boosted ADM decomposition. It 
would be also very interesting to apply the partially constrained vielbein formulation to the 
new quasi dilaton extension of massive gravity or similarly to bigravity and investigate their 
cosmological implications. 
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